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SUMMARY 

An experimental and theoretical Investigation of vapor cavitation Is made 
for a submerged journal bearing under dynamically loaded conditions. The 
observance of vapor cavitation In the laboratory was made possible using high- 
speed photography. Vapor cavitation was found to occur when the tensile stress 
applied to the oil exceeded the tensile strength of the oil or the binding of 
the oil to the surface. The theoretical solution to the Reynolds equation Is 
determined numerically using a moving boundary algorithm. This algorithm con- 
serves mass throughout the computational domain Including the region of cavlta 
tion and Its boundaries. An ADI (Alternating Direction Implicit) method Is 
used to effect the time march. A study of a rotor undergoing circular whirl 
was conducted. Predicted cavitation behavior was analyzed using three- 
dimensional computer graphic movies. The formation, growth, and collapse of 
the bubble In response to the dynamic conditions Is shown. For the same con- 
ditions of dynamic loading, the cavitation bubble was studied In the laboratory 
using high-speed photography. Movies comparing experiment with theory will be 
presented. 


INTRODUCTION 

The most common cavitation boundary condition for steadily loaded bearings 
Is the zero-pressure-gradient condition of Swift and Stieber (refs. 1 and 2). 
Ihls condition results from the assumption of continuity of mass flow at the 
boundary between the cavitation zone and the 1 1quid- f 1 1 led regions under steady 
state conditions. In many bearings, however, dynamic loads cause changes In 
ttie local film thickness, and the problem of determining the appropriate bound 
ary condition requires a condition other than zero pressure gradient at the 
boundary. The studies of Olsson (refs. 3 and 4) suggest that the usual steady 
state condition, p = dp/dx ^ 0, Is adequate for dynamic situations if the cavi 
tatlon boundary moves at a speed that Is less than half the journal surface 
speed. However, If the cavitation boundary moves at a speed that is greater 
than half the journal surface speed, a more realistic formulation of the bound 
ary condition should be developed. 
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Jakobsson-Floberg (ref. 5) and Olsson (ref. 3) formulated boundary condi- 
tions for a moving boundary that conserved mass within the cavltated region as 
well as at the boundary (commonly referred to as the JFO cavitation theory). 

It was assumed that liquid was convected through the cavltated region In the 
form of strlatlons extending to both surfaces In the film gap. The striated 
flow In the JFO theory Is necessarily a Couette flow because of a constant 
pressure assumption within the cavltated region. Jakobsson and Floberg's expe- 
rlmen al findings (ref. 5) support the constant pressure assumption. The JFO 
th‘iiy perhaps represents one of the best accounts of a dynamical theory to 
dale for moderately to heavily loaded Journal bearings and/or dampers. It Is 
an Improvement over the Swift-Stleber conditions, even for steady- state solu- 
tions, because It provides for film reformation or flllback as well as film 
rupture. Both rupture and flllback require a knowledge of the pressure gradi- 
ent and fractional film content at the Interface to determine Its location. 
Unlike a rupture boundary, the flllback boundary Is subjected to a pressure 
flow (I.e., nonzero pressure gradient). Further, the fractional film content 
at the boundary Is determined by a residual fluid within the cavltated region 
that has been released at a rupture boundary (earlier 1n time) and governed by 
the fluid transport law. Which condition (flllback or rupture) prevails, 
depends upon the relative motion of the boundary with the motion of the con- 
vected fluid normal to the boundary. Typically, under steady-state cavltating 
conditions, the rupture boundary occurs along the upstream boundary and the 
flllback boundary occurs along the downstream end. Under dynamic loading and 
nonstationary cavitation, rupture or flllback conditions may be required at 
either the upstream or the downstream boundary. To effect these conditions In 
a computational algorithm, the programming task Is exceedingly tedious and 
this discourages their Implementation In current practice. 

tlrod and Adams (ref. 6) Introduced a computational scheme that mimics the 
JFO theory. It avoids the complex programming required to trace the moving 
boundary between grid points and evaluating the required pressure derivative. 
Later, Elrod (ref. 7) modified It and presented It In more detail. Ihe algo- 
rithm Inorporates a switching function that “switches out" any pressure flow 
within the cavltated region. This automatically Introduces a cavity conforming 
to the requirements of mass continuity and the JFO theory. Elrod has compared 
results generated by the algorithm with experimental results of l.undholm 
(ref. 8) for steady-state operation of a circumferentially fed Journal bearing. 
The results (I.e., cavitation extent, load capacity, and attitude-angle versus 
eccentricity) were quite good. For steady-state applications, Lebeck et a1. 
(ref. 9) and Miranda (ref. 10) have used the Elrod algorithm. Brewe (ref. 11) 
was one of the first to apply the algorithm to a submerged journal bearing 
operating under dynamically loaded conditions. The solution to the Reynolds 
equation Is determined numerically using a control volume method. This method 
conserves mass throughout the computational domain Including the liquid-vapor 
Interface which may or may not be In motion relative to the minimum film line. 
An ADI (Alternating Direction Implicit) method Is used to effect the time 
march . 


To the authors knowledge, Jacobson and Hamrock (ref. 12) was the first 
paper to experimentally Illustrate the development of vapor cavitation In a 
dynamically loaded Journal bearing. A high-speed motion picture camera was 
used to capture the formation and collapse of the vapor bubble In the journal 
bearing . 
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The theoretical formulation presented by Brewe (ref. 11) will be used to 
describe the formation, growth, and collapse of a vapor cavitation bubble In 
response to conditions In a dynamically loaded Journal bearing. For the same 
conditions of dynamic loading, the cavitation bubble will be studied using the 
experimental apparatus of Jacobson and Hamrock (ref. 13). The development of 
the vapor cavitation will ^>e shown using the high-speed photography and these 
results will be directly c >mpared with the computer simulated results. 
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SYMBOLS 

shaft diameter, m 
eccentricity, m 

switch function (cavitation Index) 

film thickness, m 

length to diameter ratio 

lineal mass flux, kg/ms 

fluid pressure, N/m^ 

ambient pressure, N/m^ 

cavitation pressure, N/m^ 

radius of shaft, m 

radial clearance, m 

time, s 

sum of the surface velocities In x-directlon, m/s 

squeeze-velocity, m/s 

coordinate along circumference, m 

Incremental spacing along circumference, m 

coordinate normal to x, z-plane, m 

axial coordinate, m 

axial Incremental spacing, m 

liquid bulk modulus, N/m^ 

angular position of minimum film, rad 


3 


E eccentricity ratio, e/AR 

Q fractional film content In cavitation zone, density ratio (p/pc) 

In full film zone 

V dynamic viscosity. Ns/m 

fluid density within cavltated zone, kg/m 

<p angular coordinate relative to minimum film line, rad 

orbital angular velocity of Journal center about fixed point relative 
to the housing center, rad/s 

angular velocity of journal about Its own center, rad/s 

Subscript s : 
d dynamic 

s static 


APPARATUS AND TEST PROCEDURE 

The experimental setup with the high-speed camera, motors, and test rig Is 
shown In figure 1. The high-speed motion picture camera that was used could 
take up to 11 000 frames/sec with ordinary 16-mm motion picture film. In the 
experiments presented In this report the framing rate was 2000 frames/sec. 

The test bearing was lubricated by gravity feed from a can placed approx- 
imately 600 mm above the bearing. The test lubricant was a nondegassed 
Dexron-II type of automatic transmission fluid with a viscosity of 0.066 N s/m 
at atmospheric conditions. The lubricant passed from the can through a line 
and entered the bearing at the top left of the polymethylmethacrylate (PMMA) 
bearing shown In figure 1. At the right In figure 1 Is the motor that gives 
the sliding velocity to the roller, and at the left Is the motor that gives 
the squeeze velocity to the bearing. The speeds of these Identical 0.5-kW 
motors can vary continuously from 0 to 250 rad/s. 

Figure 2 shows the PMMA bearing, which can vibrate In a circular motion 
keeping the axis of vibration parallel to the axis of rotation. The special 
mechanism manufactured to keep the centerline of the hearing parallel to the 
centerline of the journal when the bearing Is vibrating Is shown In figure 2. 

It consists of three plates each coupled together with four thin plates made 
from spring steel. The spring steel plates could easily be bent In the direc- 
tion perpendicular to the plates but were very stiff In the horizontal planes 
of the plates. Because of the symmetry of the mechanism the thick plates were 
always parallel to each other when the thin plates were bent as long as they 
did not buckle. 
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The dynamic eccentricity mechanism shown assembled In figure 2 enables one 
to vary the dynamic eccentricity from 0 to 0.5 mm. Since the radial clearance 
Is 0.5 mm, the ratio of the dynamic eccentricity to the radial clearance can be 
varied between 0 and 1. The amplitude of the vibration of the squeeze motion 
Is given by the dynamic eccentricity ratio. By setting the dynamic eccentric- 
ity equal to zero and moving the journal relative to the centerline of the 
bearing, the static eccentricity can be varied between 0 and 0.5 mm, meaning 
that the static eccentricity ratio can be varied between 0 and 1. Therefore 
the offset of the center of vibration from the center of the journal gives a 
static load, and the amplitude of the vibration gives a variable rotating 
dynamic load superimposed on the static load. Both the location of the center 
of vibration and the amplitude of vibration can vary. 

Figure 3 shows four positions of the journal relative to the bearing for 
a dynamically loaded journal bearing. The static eccentricity ratio cj Is 
set at 0.34, and the dynamic eccentricity ratio «(j Is set at 0.66. The 
direction of the dynamic eccentricity rotates 90® counterclockwise going from 
(a) to (b), (b) to (c), and (c) to (d). In figure 3 the squeeze and separation 
component of motion Is clearly visible. 


EXPERIMENTAL RESULTS 

Figure 4 presents several frames from the high-speed films showing the 
formation and collapse of vapor cavitation In a journal bearing. The actual 
film provides a better Illustration; however, the figures highlight the film. 
The fixed conditions for figure 4 are L/0 = 1/4, cj = 0.340, = 0.660, 

0.32 £ 1.0, fa>s = 19.5 rad/s, wj * 92.7 rad/s, and PTFE journal surface. 

The framing rate of the camera was 2000 frames/s. The time from the first 
frame to the last In figure 4 was less than 25 ms. Indicating how rapidly the 
vapor cavitation forms and collapses. 

The sample frames filmed (fig. 4) Illustrate that In dynamically loaded 
journal bearings cavitation bubbles containing mainly oil vapor are present. 
Vapor cavitation occurred when the tensile stress applied to the oil exceeded 
the tensile strength of the oil or the binding of the oil to the surface. To 
the author's knowledge, other than our Initial report (ref. 3), this was never 
shown before. Vapor cavitation Is quite different from gas cavitation. The 
gas cavitation bubbles can be transported through the high-pressure region of 
the bearing as demonstrated In Jacobson and Hamrock (ref. 13). On the other 
hand, vapor cavitation bubbles dissolve and disappear In less than a hundredth 
of a cycle of vibration. Without a high-speed camera one would not be able to 
observe vapor cavitation. Furthermore, this type of cavitation can cause 
erosive wear (ref. 14) of adjacent surfaces, so Its presence In dynamically 
loaded journal bearings Is an Important observation. 


THEORETICAL MODEL 


Elrod and Adams (ref. 6) Implemented a moving boundary scheme that avoids 
Interface complications. The following reveals some of the underlying factors 
that led up to the algorithm which was later modified by Elrod (ref. 7) and 
used here. The conservation of mass can be written In the following way: 
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m = 
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( 1 ) 


where m represents the lineal mass flux and 1s given by 


m 




(2) 


Substitution of equation (2) Into (1) leads to the Reynold's lubrication equa- 
tion. This equation has been made applicable to the cavitation region as well 
as the full film region by Incorporating a switch function to automatically 
satisfy the boundary conditions at a moving Interface. Further, the fractional 
film content, e, has been made the dependent variable. This required giving 
0 a dual Interpretation. That Is, In the full film region, e represents the 
mass content of the film that exceeds the content that would exist If the pres- 
sure were at cavitation pressure, Pj.. In other words, 

e - p/pc (3) 

where pc 1s the density of the liquid at the pressure Pc- '"urther, p and 
e are related through the equation for the liquid bulk modulus according to 


p 


at 

3p 


m 




(4) 


or 


p = Pc + B In e 0 >1 .0 (5) 

In the cavitation region (0<1.O) determines the mass content (pc^h) which can 
manifest Itself In the form of a smeared mass or striated flow extending to 
both surfaces In the film gap. The fluid transported through the cavitation 
region In the form of an adhered film can be neglected for heavily loaded 
conditions . 


Having linked the solutions of the full film region with the solutions In 
the cavltated region via a single dependent variable, 0, a universal differ- 
ential equation Is made possible. However, a cavitation Index, or switch func- 
tion, g, was Introduced so that the resulting POE would be consistent with toe 
uniform pressure assumption within the cavltated region. The switch function 
Is defined from a knowledge of 0. Thus 


g s 0; 0 1.0 

g = 1 ; 0 1.0 


( 6 ) 


and g Is made a factor of the pressure gradient term In equation (2) so that 
the flow Is strictly Couette In the cavltated region. Expressing the lineal 

• 

mass flux, iri, In terms of 0 , and g and substituting Into equation (1) 
results In the universal differential equation obtained In Elrod and Adams 
(ref. 6); I.e., 
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( 7 ) 


^ . ^(eh) - I • ^ Bg(e)^e 
at 2 '‘w 


at 

In the fun film region the solutions of e together with eQujJJ'o" (5) 
detem^nr'h. pressures. In the cevUetlor repten. g . 0, and eguatloh (7) 

becomes 


^ . ^(eh) - 0 


( 8 ) 


governing the transport of the fluid through the cavitation region. 

Making use of Information from Brewe (ref. 11) the resulting convective 
and pressure flow equations can be expressed as 


^^""x^conv 


^^"'x^press 
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Tn this Study an alternating direction Implicit (ADI) scheme was applied. 
nfng'tSe^^rSe'rfncrrre JoSnS'SrpeJjJrm^ng a 

oerlodlc matrix using a maximum-pivot strategy to reduce the error. In the 
^ nnH half timp stpD differencing along the axial direction leads to a system 
° aM^nc il^t be exDrSssed In pure tridlagonal form. These solutions 
are^readlly obtained by using a nonpivoted Gaussian elimination procedure (tr1- 
aiLonti solver? Since the Journal and housing are aligned, the pressures 
:::?Te sy^rtric abi^tJe^ilal center Consequently, the calculations are 
only made over half the axial length of the housing. 

ThP switch function was updated after each half-time step. Occasionally. 

T.Zllrl ^hli^S^e ::iuS ruSc;irwIrn:i^a:Jo-^odSli:g%:e 7 L:nt°rj 

S££?~'Ln^.Tu‘hc^i: ^ 

was reiterated. 


7 


THEORETICAL RESULTS 


A numerical scheme (ref. 11) based upon the Elrod algorithm was devised 
to study vapor cavitation In a submerged Journal bearing for both steady and 
dynamic loading. The method conserves mass throughout the computational domain 
Including the region of cavitation. In this Investigation, the algorithm was 
checked against experimental data for a Journal whirling In a circular path 
through one orbit. 

Figure 5 Illustrates the bearing configuration along with the predicted 
pressure distribution that resulted from Brewe's (ref. 11) numerical calcula- 
tions using the Elrod algorithm. The Journal Is rotating In a counter clock- 
wise direction. Thus, *:he entrainment velocity through the minimum film Is In 
a counter clockwise direction, and by convention considered positive. The line 
of centers extends through the position of minimum film thickness. The posi- 
tive pressures seen In t'le three-dimensional plot of figure 5 are generated In 
the converging clearance; region adjacent to this line. In the direction of 
rotation, the pressures become S'lbambicnt In the diverging clearance region. 

The upstream and downstream meniscus of the cavitation boundary Is outlined 
there. The corresponding extent of cavitation Is shown by the dark shaded 
area In the three-dimensional plot. Figure 5 Illustrates the similarities and 
differences of the pressure distribution as they occur throughout the entire 
orbit. Frame (a) represents the position of the journal within the housing 
and the associated pressure distribution at the Initial Instant In time. The 
pressure buildup due to the combined squeezing and sliding motion of the Jour- 
nal Is shown for the first half cycle In frames (a) through (d). It should be 
noted In comparing frames (c) and (d), that the maximum peak pressure does not 
occur when the eccentricity Is a maximum as one might first expect. This can 
be explained by referring to figure 6 and noting that the squeeze velocity 
approaches zero and the sliding velocity approaches a minimum at the end of 
the half cycle. This offsets the effectiveness of generating Increasing pres- 
sures due to a decreasing film. The onset of cavitation o.curs at frames (b) 

25 ms Into the cycle. The squeeze component Is known to suppress the condi- 
tions for cavitation (ref. 16). Proceeding from (d) to (e), shows the Journal 
(near the minimum film line) Is separating from the bearing. The Initial 
stages of separation creates a suction effect, causing the pressure hump to 
dissipate and the vapor bubble to expand. Note that frames (b) and (e) were 
both generated at the same value of eccentricity. The essential difference Is 
that the squeeze velocity has a different sign. That Is, when the squeeze 
velocity Is negative (I.e., the clearance Is decreasing), the motion of the 
journal Is producing a squeeze effect (frame b). In frame (e), the squeeze 
velocity Is positive (I.e, the clearance Is Increasing) and produces a suction. 
Here the vapor bubble has actually crossed the minimum film line and been drawn 
Into the converging clearance space of the bearing. As the Journal continues 
to pull away from the bearing, the Increased clearance means that a much 
greater Polseullle side flow Is present to cause the collapse of the bubble. 
Ihls Is because the Polseullle side flow Is proportional to h^. in frame 
(f), the bubble has drifted downstream from the minimum film line and Is 
collapsing. 


COMPARISON OF EXPERIMENTAL AND THEORETICAL RESULTS 

The formation, growth, and collapse of the bubble 1n response to dynamic 
conditions Is compared experimentally and theoretically. High-speed photo- 
graphy revealed In figure 4 that for the conditions listed previously, the 
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vapor bubble remained visible fir a period of at least 23 ms. This represented 
approximately one-third of the full dynamic cycle of the bearing. 

The predicted life of the vapor bubble shown In figure 5 from the computer 

code was 32 ms. A part of this difference arises because this measure of 

bubble life Inherently contains a certain duration of time for which the com- 
puter Indicates cavitation although In reality It would be Invisible to the 
eye. This Is probably not the major contributing factor to the discrepancy, 
however. The pressures calculated for these conditions were of the order of 
10^ N/m^. A plastic (polymethylmethacrylate, PHMA) was used as a housing mate- 
rial for viewing purposes. It Is thus most likely that deformation and temper- 
ature effects In the experiments are responsible for the discrepancy, not 
accounted for by the numerical method. Further, the experiment allowed the 
eccentricity to reach a maximum value of 1.0, whereas a certain amount of 
numerical Instability was observed when the eccentricity exceeded 0.98. Con- 
sequently, the Input conditions for the computations were not Identical to the 
Input conditions of the experiment. 


CONCLUSIONS 

A theoretical and experimental Investigation Is made of the evolution of 
a vapor-bubble for a submerged Journal bearing under dynamically loaded condi- 
tions. The theore!.1cal solution to the Reynolds equation Is determined numer- 
ically using a control volume method (Elrod algorithm). This method conserves 
mass throughout the computational domain Including the liquid-vapor Interface 
which may or may not be In motion relative to the minimum film line. An ADI 
method Is used to effect the time march. 

The occurrence of vapor cavitation was clearly evident In both the theo- 
retical as well as the experimental results presented. Vapor cavitation occur- 
red when the tensile stress applied to the oil exceeded the tensile strength of 
the oil or the binding of the oil to the surface. The physical situation caus- 
ing the vapor cavitation was a squeezing and sliding motion within a bearing. 
Good agreement was found between theory and experiment. A prediction of 32 ms 
bubble life for nonstationary cavitation was obtained from the dynamic theory. 
Using high-speed photography for the same set of conditions, a bubble life of 
23 ms was observed In the laboratory. This discrepancy was attributed to 
deformation and temperature effects as well as some uncertainty In correlating 
actual life with duration of visibility. Furthermore, the maximum eccentricity 
In the experiments was 1.0. Numerical Instabilities prohibited eccentricities 
larger than 0.98. 
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Figure 1. - Test apparatus. 
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Figure 2. - PM^tA'tube and eccentric. 
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16 Abstract 

An experimental and theoretical Investigation of vapor cavitation Is made for a i 
submerged Journal bearing under dynamically loaded conditions. The observance | 
of vapor cavitation In the laboratory was made possible using high-speed photo 


graphy. Vapor cavitation was found to occur when the tensile stress applied to 
the oil exceeded the tensile strength of the oil or the binding of the oil to 
the surface. The theoretical solution to the Reynolds equation 1s determined 
I numerically using a moving boundary algorithm. This algorithm conserves mass 
I throughout the computational domain Including the region of cavitation and Its 
boundaries. An ADI (Alternating Direction Implicit) method Is used to effect the 
time march. A study of a rotor undergoing circular whirl was conducted. 

Predicted cavitation behavior was analyzed using three-dimensional computer 
graphic movies. The formation, growth, and collapse of the bubble In response 
to the dynamic conditions Is shown. For the same conditions of dynamic loading, 
the cavitation bubble was studied In the laboratory using high-speed photography. 
Movies comparing experiment with theory will be presented. 
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